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Abstract. Consider the non-compact billiard in the first quandrant bounded by 
the positive s-semiaxis, the positive y-semiaxis and the graph of f(x) = (x + l)~ a , 
~ff-)~ ■ a E (1, 2]. Although the Schnirelman Theorem holds, the quantum average of the 

position x is finite on any eigenstate, while classical ergodicity entails that the 
■ classical time average of x is unbounded. 
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The purpose of this note is to exhibit a simple example of a chaotic system in which 
the long time classical and quantum evolutions can be proved to be qualitatively dif- 
ferent even though the Schnirelman Theorem (sometimes called quantum ergodicity: 



see jSc] and also [Hc!V| . jDTTTj . jHMRj . [Zi], [ZZ]) holds. 

Let Q be the planar domain bounded by the positive x-semiaxis, the positive y- 
semiaxis and the graph of f(x) = (x + l)~ a , a £ (1, 2]. Notice that Area(Q) < +oo. 
Let 0* : QxS 1 — ► Q x S 1 be the dynamical flow corresponding to the billiard motion 
in Q. This is the Hamiltonian flow of a particle of energy 1 that moves freely in 
Q and performs totally elastic collisions at the boundary. By |Le2j 0* is ergodic 
w.r.t. the normalized Liouville measure dv := (2tt Area(Q)) -1 dxdydO. (Ergodicity 
there is proved for some Poincare sections and then extended to the billiard flow; see 
also |Le3j for a general class of ergodic billiards of the like, of finite or infinite area, 
including some that qualify as further examples for the present paper.) Moreover, 
as we show in the Appendix below, the Lyapunov exponent of the flow is positive. 
Therefore this system is chaotic according to general consensus. 

In quantum mechanics, the corresponding system is defined by the Schrodinger 
operator 

H:=-h 2 A D , D{H) = H 2 {Q)nHl{Q) C L 2 (Q), (1.1) 

where Ad is the Laplacian with Dirichlet boundary conditions at dQ. It is an old 
result of Rellich that Ajj is self-adjoint with compact resolvent. Hence Spec(if) is 
discrete. (More recent results on Spec(A£>) include |Taj . |Roj . [51] . |vdB| . |Daj ) . We 
denote Ej(h) = h 2 £j, j £ N, the eigenvalue of H corresponding to the normalized 



SANDRO GRAFFI MARCO LENCI 



eigenfunction if)j. The order is such that j i— > Ej is non-decreasing and {V'jjjeN is a 
complete orthonormal system. 

The classical time average of the coordinate x is 



X(x, y,9):= lim ^ T X o <f>\x, y, 9) dt, (1.2) 

where X(x,y,8) := x. In principle one does not know whether this limit exists for 
almost all (x,y,8), as X ^ L l (Q x S* 1 ) and Birkhoff's Theorem does not apply. 
However, Proposition 11.11 below will override this question. 

Tolerating a little abuse of notation, the quantum expectation of the position 
operator X on any u G L 2 {Q) is instead 

(X)(u) := (u, Xu) = / x\u(x, y)\ 2 dxdy, (1.3) 



JQ 

The comparison result we want to prove can be formulated as follows. 
Proposition 1.1. 

(a) The classical average is infinite, i.e., for v-a.e. initial condition (x,y,6): 

X(x,y,6)=+oo. (1.4) 

(b) The eig en functions of Ad are super-exponentially localized in the following 
sense: Denote 

^(x) := / \^(x,y)\ 2 dy, (1.5) 

the probability density of the position in the x-direction. Then one has: 

V 7 >0, £,■(*) = (e~n (1.6) 

as x — > +oo. Therefore, in particular, the quantum expection of X on every 
eigenstate is finite: 

r r+oo 

(X)(ipj) :— / x\ijjj(x, y)\ 2 dxdy = / x£j(x)dx < +oo. (1.7) 
Jq Jo 

(c) The Schnirelman Theorem holds. As a corollary, there is a density-1 sequence 
{jn}, i-e-, a sequence with the property 

lim = 1, (1.8) 

k— >+oo K 



such that 



lim (X)(^-J = +oo. (1.9) 

n— »+oo 
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Remarks 

(1) Proposition !!. 11 (b) implies that (F(X))(u) < +00 for any (normalized) state 

00 

u := ^^ajipj with j 1 — > \dj\ decaying sufficiently fast. Here F G L} oc (M) is 
i=i 

any x-dependent observable exponentially bounded as x — > +00. In fact, a 
Cauchy-Schwartz inequality shows that 

\(u,F(X)u)\ < \\F(X)u\\ (1.10) 

00 / >"+oo \ 1/2 

^ 5Zl a ^(^y |-F(x)| 2 ^(x)rfxJ < +00. 

If we replace F(X) by the corresponding Heisenberg observable F(X)(t) := 
e iHt/h F(X)e- iHt/h , clearly the bound (fOTljl holds uniformly in i: 

\(u,F(X)(t)u}\ = \(e~ iHt/h u,F(X)e- iHt/h u}\ (1.11) 

1/2 



This last estimate shows that the quantum evolution is also super-exponen- 
tially localized. 

(2) Therefore, for the same u as above, one also obtains convergence of the time 
average of the Heisenberg observable X(t): 



POO) 



lim - / (u,e iHt/h Xe- im/h u) 

a*a k {ijj j: Xip k ) < +00. (1.12) 



i, kef 



This Heisenberg time average (X) is the quantity to be directly compared 
with the classical time average X. 
(3) Assertion (b) shows that quantum mechanics localizes the unbounded clas- 
sical chaotic motion as soon as h > 0. This phenomenon is an example of 
quantum suppression of classical chaos. Its physical origin is clear: when 
the quantum particle visits the deepest recesses of the cusp, its ?/-position is 
very well determined. Therefore, by the Uncertainty Principle, its y-velocity 
must be extremely spread out, and this cannot occur at finite energies. (This 
crucially depends on the Dirichlet boundary conditions: for the Neumann 
Laplacian the situation is quite different |DSj .) One might think of this phe- 
nomenon as the opposite of the tunnel effect: the classical particle is much 
more likely to "penetrate" the cusp than the quantum one. 



4 



SANDRO GRAFFI 



MARCO LENCI 



(4) The classical limit is naturally defined as the joint limit j — > +00, h — > 0, such 
that -Ej(^) = 1 (1 being the fixed value of the energy). Since the quantities at 
hand here do not depend on H, the second limit can be forgotten. Obviously, 
(jl.6j) - (jl.7j) are not uniform as j — > +00. 

(5) In part (c) the Schnirelman Theorem is formulated for the first time for 
systems whose classical energy surface is not compact (although one does 
not really need more than |ZZj to prove it in our case — see below). As is 
well known, the Schnirelman Theorem is a purely asymptotic statement. 
The present example shows moreover that this asymptotic statement cannot 
exclude that quantization can turn a classical behavior at infinity into a 
behavior of a completely different nature. Hence, in general, it might appear 
physically misleading to call quantum ergodicity the simple validity of the 
Schnirelman Theorem. 



the assertion follows follows directly from classical ergodicity |Le2j and the following 
easy lemma. 

Lemma 2.1. Let (V, v) be a probability space and </>* a flow on V that preserves the 
measure v . If the measurable function g is bounded below and (V,^, v) is ergodic 
then, for v-a.e. z G V , 



whether g is integrable or not. 

Proof. If g is integrable there is nothing to prove by ergodicity. Otherwise J v g dv = 
+00, because g is bounded below. In this case, for m 6 N, set g m := i m o g, with 
i m {x) '■= x for x < m and i m {x) '■— m for x > m. By the Birkhoff theorem, there is 
a full-measure set B m over which the time average of g m exists and is equal to its 



2. Proof of the Proposition 



(a) Since clearly 




(2.1) 




(2.2) 



LOCALIZATION IN INFINITE BILLIARDS 



phase average. Therefore, for z G B := C\ B m (still a full-measure set), 



lim inf — 



mgN 

f z) dt, > lim inf — 

fl'm diA 



g m o ( f) t (z) dt 



(2.3) 



By monotonic convergence, the sup of the r.h.s. in m is +00, while the l.h.s. does 
not depend on m. This proves at once that the limit in T exists and is +00. □ 

(b) For the proof of this part, let us drop the subscript j from all notation. Hence 
ip is continuous on Q, infinitely smooth in its interior, and such that 



-(d 2 x + d 2 U(x,y) = £i;(x,y), 



(2.4) 



with ip\dQ = 0. Without loss of generality if) is real. Recalling definition (jl.5j) one 
easily checks, by repeated differentiation inside the integral, that 



/•/(«) rf( x ) 
e(x) = / dl^dy > 2 / 
Jo Jo 



ipd^ipdy. 



(2.5) 



Now, think of —d 2 as the 1— dimensional Laplacian on [0, a] with Dirichlet boundary 
conditions. Then, in the sense of the quadratic forms, — d 2 > (71 /a) 2 . Therefore, 
multiplying ()2.4j) by tp, integrating over y, and using that lower bound, we obtain 



/(*) 



ipdlipdy 



Plugging into (j2^ 



e(x) > 2 



< 



7T 



ip 2 dy + 
2 



/(*) 



7T 



fix) 



ip 2 dy. 



(2.6) 



(2.7) 



for any 7 > 0, provided x is large enough depending on 7. This means that either 
£(x) > Ce lx or ^(x) < Ce~ rx . Since t/> is an eigenfunction, the first possibility 
cannot occur. (jl.6|) is thus proved. 

(c) The version of the Schnirelman Theorem that can be verified in our case is the 
following: For any pseudodifferential operator A of order 0, compactly supported in 
the position variables, there exists a density-1 sequence {j n } such that 



lim (if>3 n ,Aif> jn ) 



a dv, 



f2.8) 



QxS 1 



where a is the principal symbol of A. The proof is a strightforward check that all 
arguments of |ZZj valid for compact billiards hold true in this case too. We limit 
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ourselves to remark that the crucial fact is that Area(Q) < +00, and in particular 
that the system is recurrent (see |Lelj and references therein for better results in 
this direction). 

To prove (jl.9j) we use an argument that is somewhat similar to Lemma 12.11 Let 
{X m } m6 N be an increasing sequence of functions on Q x S 1 such that: X m is smooth 
and depends only on x, supp(X m ) C {x < m}, X m < X and, pointwise, X m — > X for 
m — > +00. With the usual abuse of notation, we denote X m also the corresponding 
multiplication operator on L 2 (Q). 

For every m there is a density-1 sequence {j n m ^} n &N = : such that ()2.8|) holds 
for A = X m (thus a = X m ), using that sequence. Now fix pi = and, for m > 2, 
consider the following recursive definition: Select a sufficiently large p m > p m -i so 
that 



Vjl m) > Pn 



(ip( m ),X m ip (m) ) - / X m dv 



xS 1 



< 3. (2- 9 ) 



m 



and 



#(q (m - 1) n[ Pm _ 1 , Pm )) ^ i_ (2 iQ) 

Pm Pm~l 

(here we have identified with its image). ()2.10|) can always be verified because 
(j( k > fl [pk, +00) has density 1. 

Now set cr := ^ ) (a (m ' ) fl [p m ,p m+ i)). Think of a as a sequence and label its 

m gN 

elements {jn}. By ()2.1()j) . <r has density 1. Moreover, by (j2.9j) . 

(ip^X^) > (i> jn ,X m ^ jn ) > [ X m dv-±-, (2.11) 



J QxS 1 m 

it jn £ [pm^m+i)- Therefore the liminf of the l.h.s. in n corresponds to the liminf 
of the r.h.s. in m: the latter is infinity by monotonic convergence. 

3. Appendix: Positivity of the Lyapunov exponent 

In this section we prove that, as claimed in the Introduction, 0* has a positive 
Lyapunov exponent A (and thus a negative exponent —A, as 0* is volume-preserving 
and the exponent in the direction of the flow is obviously 0). 

For z := (2, y, 6) G Q x S 1 , denote E c z C T 2 (Qx S 1 ) the one-dimensional subspace 
of the tangent space at z generated by dt4> t (z), i.e., the direction of the flow in z. (The 
inconsequential complications arising when (x,y) G dQ are ignored). It is known 
that D(f) 1 preserves the splitting E c © (-E C )" L (see [SC a — this corresponds to the fact 
that the wave- front of a light wave is always orthogonal to its rays). Therefore the 
stable and unstable directions at z (denoted E s z and E™, respectively) are to be 
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looked for within (Ef) L . To verify their existence we use the results of |Le2j about 
A4, the Poincare section given by all unit vectors (line elements, in the language of 
[S]) that are based in dQ and point towards the interior of Q. We know that T, the 
Poincare map of A4, has stable and unstable directions, denoted F s ^: we will use 
these to define E s ^ everywhere (more precisely, in a full-measure subset of Q x S 1 ). 

Consider a non-singular element z G A4, that is, a line element whose future 
orbit never hits dQ tangentially or intersects the vertex at (0,1). (Singular line 
elements can be neglected because they form a measure-zero set in A4, and thus in 
Q x S 1 , relatively to their respective invariant measures). Since z is non-singular, 
E c z <f_ T z Ai and so IT 2 : T Z A4 — > (-E^) -1 , the projection in the direction of E°, is 
non-degenerate. We define Ez := H z Fz . Finally, for t > 0, DftEl^ define the 
stable and unstable directions at all points in the future orbit of z. It is easy to see 
that this definition is unambiguous. 

For a non-singular z (not necessarily in A4), take v u G E%. The sought Lyapunov 
exponent is 



1 I Drh T (v u )\ 1 
A(*):= lim -log 1 )l = lim -log DfA^ . (3.1) 



Let us introduce the following notation: r(z) := min {i > | (^(z) G dQ x S 1 } is 
the (forward) free path of the line element z, t~(z) is the backward free path, 
analogously defined, and the closest line element in the 

past orbit of z that is based in the boundary. Define 









E u . 

z' 



(3.2) 



where D<ft t+ means lim D<j) t+£ . (This is needed because 0* is discontinuous at 
dQ x S 1 .) By construction, g is constant along the rectilinear parts of a given 
orbit. Furthermore, due to the composition properties of the differential, 



1 



T 



X(z)= lim - / go4>\z)dt. (3.3) 

In fact, in (|3.1|) and ()3.3|) . the two functions that are averaged in T can only differ 
in the first and last segment of of (z)}J =Q . This difference gets washed away as 
T -> +oo. 

We claim that g > 0. In order to see this, let us go back to the Poincare map T. 
For z G Ai, we have seen that 

D^ z)+ \ {Ec ^=TiT {z )DT z Yi-\ (3.4) 



8 SANDRO GRAFFI MARCO LENCI 

But it is a fact that DT\ pu is strictly expanding w.r.t. to a certain metric in TA4. 
This metric is called in |Le2j the increasing metric and is defined as |f|i nc := |IIt>|. 
This fact and (|H.4|) show that D<p T+ \ Eu is strictly expanding which, in view of (JS2J), 
proves the claim. 

Finally we apply Lemma f2. II to conclude that, for u-a.e. z G Q x S 1 , 

\(z) = I gdv>0. (3.5) 

JqxS 1 
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